Transient dynamics of Light propagation in EIT medium and hidden symmetry of 

multi-bit quantum memory 
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We investigate the transient phenomenon or property of the propagation of an optical probe field 
in a medium consisting of many A-type three-level atoms coupled to this probe field and an classical 
driven field. We observe a hidden symmetry and obtain an exact solution for this light propagation 
problem by means of the spectral generating method. This solution enlightens us to propose a 
practical protocol implementing the quantum memory robust for quantum decoherence in a crystal. 
As an transient dynamic process this solution also manifests an exotic result that a wave-packet of 
light will split into three packets propagating at different group velocities. It is argued that "super- 
luminal group velocity" and "sub-luminal group velocity" can be observed simultaneously in the 
same system. This interesting phenomenon is expected to be demonstrated experimentally. 
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I. INTRODUCTION 

The recent experiments [1,2] have demonstrated many 
exotic natures of light pulse propagation in the solid 
and gas systems with electromagnetically induced trans- 
parency (EIT) [3], such as the phenomena of "super- 
luminal group velocity" and "sub-luminal group veloc- 
ity". Physically they reflect the quantum coherence ef- 
fects to correlate the quantum fluctuations [4]. The 
concepts of atomic coherence and interference have also 
been applied to lasing without inversion [5] and the en- 
hancement of linear or nonlinear susceptibilities [6,7]. At 
present the sub-luminal phenomena including the light 
stopping in an atomic medium have been extensively 
studied from various points of view [8]. The recent ex- 
periments have also shown the ultra-slow group velocity 
of light in the solid state systems, such as 3-mm thick 
crystal of Y 2 SiO^ [9]. Most recently, it is found that 
when a few rubidium atoms are loaded into high-Q opti- 
cal micro-cavity, the " super-luminai" and "sub-luminal" 
phenomena can be observed in a refined version of the 
above mentioned experiments [10]. 

With possible applications in quantum information 
field developed in the past ten years [11], some studies on 
the sub-luminal problem are closely associated with an 
ideal and reversible transfer technique for the quantum 
state between light and metastable collective states of 
matter [12]. Theoretically, the basic idea is based on the 
control of light propagation in a coherently driven 3-level 
atomic medium. The exciting fact that the group velocity 
is adiabatically reduced to zero [12] means the possibil- 
ity of proposing a more practical protocol to store and 
transfer the quantum information of photons in the col- 
lective excitations in an ensemble of A-type 3-level atoms 
[13,14]. This may open up an interesting prospective for 
quantum information processing. The studies show that 
the excitations of matter coupled by light are more sta- 



ble under some circumstances and thus form the so-called 
" dark-state polaritons" . Physically speaking, this is the 
essence of the problem. Extended to the multi-atom case, 
the above idea about adiabatic transfer of the state of a 
single photon into that of an individual atom [15] pro- 
vides a conceptually simplest approach for the implemen- 
tation of the quantum memory of photon information. 
Technically, this approach combines the enhancement of 
the absorption cross section in multi-atom systems with 
dissipation-free adiabatic passage techniques. In addi- 
tion, this significant investigation motivates the protocol 
for long distance quantum communication [21] based on 
atomic ensemble. 

Most recently we have deeply investigated the robust- 
ness of this kind of quantum memory [22]. To avoid 
the spatial-motion induced decoherence in an ensemble 
of free atoms, we have naturally proposed a protocol that 
each A-type atom is fixed on a lattice site of a crystal [17] . 
As quantum memories robust for the spatial-motion in- 
duced decoherence, the quasi-spin wave collective excita- 
tion of many A-type atoms forms a two-mode exciton sys- 
tem with a dynamic symmetry (or a hidden symmetry) 
depicted by the semi-direct product algebra SU(2)®ti2 
(/12 is an algebra of two mode boson operator). Physi- 
cally, this hidden symmetry guarantees the stable spec- 
tral structure of such dressed two-mode exciton system 
while its decoherence implies a symmetry breaking. If 
one only considers quantum memory, one can focus on 
the case of single mode light field. However, if one con- 
siders the propagation of light pulse in this EIT medium, 
it is then necessary to concern the multi-mode light field 
since a light pulse can be understood as a superposition 
of the components of different frequencies. The aim of 
the present and subsequent papers is to analyze the tran- 
sient phenomenon or property of the light propagation in 
this EIT medium, especially to emphasize the role of the 
generalized hidden symmetry. 
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In this paper, our proposed system still consists of the 
quasi-spin wave collective excitation of many A-type 3- 
level subsystems as in ref. [17] . The meta-stable state still 
interacts with an exactly resonant classical field, but the 
approximately resonant quantized light field, which cou- 
ples to the transition between the excited state and the 
relative ground state, is no longer of single mode. We will 
prove that, in the large N limit, one needs to introduce 
a pair of exciton operators for each mode of quantized 
light field. These realize the infinite boson algebra local 
in the mode space (or the frequency domain), but the 
collective quasi-spin operators intertwining between the 
meta-stable states and the excited ones generate a sin- 
gle global £{7(2) algebra. With the help of this hidden 
symmetry and the corresponding spectral generating al- 
gebra method, we exactly obtain the dressed spectra of 
the total system formed by the infinite two-mode exci- 
tons coupling to a multi-mode quantized electromagnetic 
field in the large N limit. The exact solutions for the 
eigen-states obtained in this way also include the multi- 
mode dark states extrapolating from photon to one mode 
quasi-spin wave exciton. Actually, these dressed states 
describe the polaritons only coupling a quasi-spin wave 
exciton to a photon for some special case. With the 
EIT mechanism the external classical field can be arti- 
ficially manipulated to change the dispersion properties 
of the quasi-spin wave excitonic medium dramatically so 
that the quantized probing light can propagate in exotic 
ways. In this way the quantum information can be co- 
herently stored and transferred among the multi-mode 
cavity photon and the multi-mode exciton system. Our 
studies in this paper are substantially related to the hid- 
den symmetry and its breaking. The hidden symmetry 
leads to an exact class of solutions for the light propa- 
gation problem, showing that a wave-packet of light will 
split into three wave-packets which propagate at different 
group velocities, namely, the " super-luminal group veloc- 
ity" , the " sub- luminal group velocity" and the usual light 
group velocity. 



II. COLLECTIVE EXCITATION OF 
THREE-LEVEL MEDIUM 

We consider the transient process for the weak light 
propagation in a medium consisting of many three-level 
subsystems. It can be an ensemble of N free A-type 
atoms, or a crystal with N lattice sites attached by N 
A-type subspaces. In recent years, the similar exciton 
system in a crystal slab with spatially fixed "two- level 
atoms" has been extensively discussed with the empha- 
sis on fluorescence process and relevant quantum deco- 
herence problem [16], [18]. 

As shown in Fig. 1, the A-type subsystem possesses 
an excited state \a), a ground state \b) and a meta-stable 
state |c). The transition frequency w QC from \a) to \c) 



of each atom is resonantly driven by a classical field of 
Rabi-frequency Q. The transition frequency w a h from |a) 
to | b) is coupled to a multi-mode quantized field with the 
annihilation operator a^, and the coupling constant gk for 
the optical mode of wave vectors k. Strictly speaking, a 
multi-mode field can't be resonantly coupled to an atomic 
transition since it is a superposition of many components 
of different frequencies. We should emphasize here that 
we only consider the case that the quantized field is a 
Gaussian wave packet in the frequency domain and the 
center frequency ujq just equals U! a b- When the frequency 
width of the wave packet is very small compared with 
<x>o (Aw <C cl>o), from the viewpoint of approximation it 
is reasonable to assume that the multi-mode quantized 
field couples resonantly to the atomic transition from |a) 
to | b). Therefore, the interaction Hamiltonian of total 
system reads 

N 



H = J2Y1 9kdk exp(ik • r^oif, 

j = l k 

N 

+Q^exp(iq- rj)o£ c + h.c, 



(1) 



where Tj (j = 1, 2, • • • , N) denotes the position of the 
j*' 1 subsystem, N the total atomic number, k the wave 
vector of quantized light of k" 1 mode and q the wave 
vector of the classical light field. The flip operators 
cr^p = \a)jj(/3\ (a, (3 = a,b, c) for a ^ (3 define the quasi- 
spin between the given levels a and /?. The coupling 

constant gk — ~P\J^^^ depends on the matrix element 
p of the electric dipole moment between \a) and \b). For 
simplicity, gk and O are considered as real without loss 
of generality. As a matter of fact, we need not require 
the resonances if we only consider the light propagation 
in such EIT resonance [24] [25]. 
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FIG. 1. Configuration of the quantum memory with A-type 
atoms, (a) attached on lattice sites of crystal, (b) resonantly 
coupled to a control classical field and a quantized probe field. 

Generalizing the definition of the excitation operators 
of single optical mode in the Ref. [17], we introduce a 
class of collective excitation operators of optical multi- 
mode 



(2) 



with respect to the transition from \b) to \a) for each wave 
vector k. Correspondingly, the collective virtual transi- 
tion from | b) to |c) can be described by another class of 
collective operators 



1 N 



N U 



(3) 



where Q(k) = k — q means the momentum conserva- 
tion in the virtual process of collective transition from 
| b) to |c). These collective operators create the general 
collective excitations 



m. 



n )=ri(^f=f4 mfc ^)ib) 



defined for the set of multi-indices 

m = (7711,7712,...), n = (771,712,...), 

and the collective ground state 

|b) = \b,b,---,b) 

with all N atoms staying in the same single particle 
ground state \b). For example, the single particle excita- 
tions \l a )k = A' k \b) and |l c )fe = C||b). These excitations 
are easy to understand. Indeed, it is obvious that, from 
the ground state |b), the first order perturbation of the 
interaction creates the so-called one exciton quasi-spin 
wave state 



N 



3-th 



,b), 



\lc)k 



1 



N 



Ye iC « k > r *\b,b, 
N f-f 



i-th 



,b). (4) 



Physically, in the large N limit under the low excita- 
tion condition, that is, there are only a few atoms oc- 
cupying the states \a) or |c) [19] and the population of 
state | b) is approximately N — > oo, the above two classes 
of quasi-spin wave excitations behave as two classes of 



bosons since they have the following bosonic commuta- 
tion relations 



[A k ,A\,] =5 k , k ', 
[Cfc, C k ,] = Sk,k', 



(5) 



To prove the above basic commutation relations, we first 
calculate 

^ A U = = 4 E e-^'-^ (ai b - <J 



1 N 

^E^ (k " k) - : 



N 

Here we have considered that 

N N 



E e - i(k '" k) - ri (^ -O - E e - i(k '- k) ' r M 6 (6) 
j=i j=i 



while 



N 



j=i 



(7) 



in the case of low excitation. For a crystal with regular 
lattice structure, according to the theory of solid state, 
we should have 



^Ee- i(k '- k) '^=^. 



(8) 



j=i 



In fact, this result can be extended to the case of generic 
medium. For a medium with a scale L much larger than 
TET ' 1PJ (this is valid for normal case) , the summation 
over Yj can be replaced by integral and it follows that 



j=l 



(9) 



since the medium volume V^ e — * oo. The other commu- 
tation relations can be proved in a similar way. 



III. PARALLELISM QUANTUM MEMORY WITH 
HIDDEN SYMMETRY 

As argued above physically the operators Ak and Ck 
depict the collective-excitation processes of bosonic type. 
The importance of these operators lies in that they define 
an invariant subspace 



3 



Vc ■ span{\m, n)|m = (m^rr^...), n = (ni ; n2,...)} (10) 



It follows that the dark-state polariton operators 



for the interaction Hamiltonian (1). This means that 
driven by this Hamiltonian, any collective state from Vc 
involves to a new collective state still in Vc- Thus we can 
use the collective excitations as basic blocks to describe 
the quantum dynamic process with a hidden symmetry. 

Let us introduce the following additional collective op- 
erators concerning the transition from |c) to \a) 



N 

£ 

j=i 



e-^ J o ,T + = (T_)t 



(11) 



which is resonantly driven by a classic light. Together 
with the third collective operator 



T 3 



EKa-<)A 
j=l 



they generate SU(2) algebra globally. The commutation 
relations between SU (2) and the collective operators A k 
and Cfc are easy to calculate. The non- vanishing commu- 
tators are as follows: 



[T-,C k 



-A k ,[T + ,A k ] = -C k , 



T+,Ct} =A{,\T_,A{ 



Cl- 



in) 



The above close commutation relations and the expres- 
sion of the interaction Hamiltonian in terms of these ex- 
citation operators: 



H = J2 9k^Na k A{ + QT+ + h.c. 



(13) 



show that there is a dynamic "group" (algebra) Gd for 
the light excited system, which is generated by A k , C k , 
A\, C\ (k = 1,2, ...), T± and T 3 . Using E to denote the 
Heisenberg-Weyl algebra generated by A k , C k , A[ and 
C, (k = 1,2, ...), we observe that the dynamic "group" 

Gd = SU(2) ® S is a semi-direct product of SU{2) and 
5 because 



[SU{2),E]cE. 



(14) 



Actually the above dynamic symmetry of Gd is a straight- 
forward generalization of the symmetry of the single 
quantized optical mode. With this symmetry the Hamil- 
tonian H can be diagonalized in an elegant way by means 
of the spectrum generating algebra method [23] . 

To this end we calculate the commutators of H with 
the generators of Gd and light field operators respec- 
tively: 



[c k ,H] = nA k ,[A k ,H] = nc k , 

[a k ,H] =g k ^A k . 



(15) 



D k = a k cos 6 k - C k sin 9 k (16) 
commute with the Hamiltonian H for the 9 k satisfying 



tan 9 k 



g k VN 
n 



(17) 



It is obvious that the dark-state polariton operators sat- 
isfy the bosonic commutation relations and define new 
dressed excitations mixing the electromagnetic field and 
collective excitations of quasi spin wave. Especially these 
new excitations are stable since 



[D k ,H}=0. 



(18) 



For the construction of the complete collective space 
dressed by the quantized light fields, another ingredient 
is the bright-state polariton operators 



B k = a k sin 6 k + C k cos 6 k 



(19) 



satisfying 



[D k ,B k 



0. 



It also extrapolates from the light field of a k to the ex- 
citon of C k when one adiabatically changes 9 k from ^ 
to zero. Evidently, the product state |0) = |b) ® |0); 
is an eigen-state of H with zero eigen- value where |0)/ 
is the vacuum of the electromagnetic field. So we can 
construct a degenerate class of zero eigen- value-states or 
dark states 



\d(n))^\d(n,t))=H^L=D{ n *\0) 



(20) 



where the adiabatic time-dependence originates from the 
change of Rabi-frcquency for the artificially-controlled 
classical field. The fact that dark states are cancelled 
by H means they can trap the electromagnetic radiation 
from the excited states. Physically this is due to quantum 
interference cancelling. 

With the help of the above hidden symmetry we can 
also construct the energy spectra with non-zero eigen- 
values by using the bright-state polariton operators. We 
start from the derived commutation relations from the 
above equation (15) 



[H,Bl] = Q k Al 
[H,Al] = G k Bl 



where 



O k = 



(21) 
(22) 

(23) 



This prompts us to consider the two commuting quasi- 
boson operators 
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Qk± = -^={A k ±B k ). 



(24) 



quantum information of photons is described by the su- 
perposition state 



As we have the commutation relations 

[H,Ql ± ]=±e k Ql ± , (25) 

the operators Q k ± can serve as ladder operators to gen- 
erate the spectra of H. In fact they produce the dressed 
states of the multi-mode exciton system 

|e(m,s,n)> = l[-^=Q^Q^\d(n)) (26) 

as the eigen-states of H with the eigen- values 

E{m,s) = ^2(m k -s k )e k , (27) 

k 

of infinite degeneracy degrees. 

Now we consider the relation between the degeneracy 
of spectrum and the hidden symmetry. We notice that 
the infinite degeneracy of |e(0,0,n)) = |d(n)) results 
from the generalized translation symmetry described by 

T fe = cxp(a k Dl -a* k D k ) 

= V ak {a k cos6 k )V Ck {-a k sm6 k ), (28) 

where 

T> a (a) = exp(aa' — a* a) 

is the usual coherent state generator. In fact, for 
the original Hamiltonian (1), the unitary transforma- 
tion T> ak (a k cos8 k ) changes the term g k \f~Na k A\ by 
gkVNA\\a k cos9 k while T>c k (— a k sin 6^) changes the 
terms flT + by — g k \J~NA\.a k sm9 k since 

T>- k (a k cos6 k )a k V ak (a k cos6» fe ) = a k + a k cos6 k , 

V^} (-a k sin 6 k )T + T>c k (-a k sin6 k ) =T + - A\a k sm.9 k . 

These two changes cancel each other. Besides this degen- 
eracy for different n, due to the relation 

[H,Q{ + Ql_}=0, (29) 

there are a large additional class of dark states of zero- 
eigen-value. They can be constructed by acting Q k+ Q k _ 
repeatedly on |d(n)). 

The above arguments suggest that, compared with the 
single mode exciton system the dressed multi-mode ex- 
citon system generated by A\,, Bj, and a\ can serve as a 
quantum memory with more advantages. It is very inter- 
esting that all the parameters 9 k for different modes of k 
can be well controlled to change from zero to infinity by 
a unique extremely slow adiabatic parameter £1 = Q(t) 
varying from a very large value to zero. As usual, the 



|L(0)) = £ c |n>, = £ c J] ^4" k 10}* (30) 

n n fc v 

where |n) ; is a multi-mode Fock state for the light field. 
In the case O(0) — ► oo, whatever the values of the pa- 
rameters are we have sin0 fc = so |d(n, 0)) = |b) ® |n) ; 
and the initial state of total system can be expressed as 
a superposition of the adiabatic dark states |d(n, 0)) at 
time t = : 

ism = |b> ® Yl c " i n ) ; = E c » i rf ( n < °)) ( 31 ) 

n 1 1 

By adiabatically changing Q to zero at time T, each com- 
ponent in the above superposition varies in this way: 

|d(n,0)) -» |d(n,T)) = |n) c ® |0>,. 

Thus one extrapolates |5(0)) to a pure multi-exciton 
state 

|5(T)) = ^c n |n) c ®|0) I) (32) 

n 

where 

|n) c = n^^|b) (33) 

is the multi mode exciton state describing the collective 
excitation. 

It is worth pointing out that, the above observation 
based on the dynamic symmetry does not depend on the 
adiabaticallity of manipulation on the external param- 
eters at all. In our example, because the dark states 
don't contain any atomic excited state \a), the sponta- 
neous emission is then forbidden even for a non-adiabatic 
manipulation. By the way we also point out that the 
quantum parallelism of the multi-mode exciton quantum 
memory gives rise to a very convenient encoding process. 
For a decimal system number n there is a unique binary 
representation: 

L-l 

n=J2 n * T ( 34 ) 

i=0 

where rii = 0,1. Then we can encode it in a special 
multi-mode Fock state 

|n); = \n a n 1 n 2 ,...,n L - 1 ) (35) 

of light field. The multi-mode Fock state |n)/ can be 
represented by an array of the many-photon qubits, with 
at most one photon in each photon qubit. It is easy 
to prepare such photon state in flying (photon) qubit 
firstly and then store the decimal system number n in the 
above multi-mode exciton quantum memory as a single 
component many-exciton state |n)c- 
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IV. TRANSIENT DYNAMICS OF LIGHT 
PROPAGATION 

In the previous sections, making full use of the hid- 
den symmetry, we not only construct two classes of dark 
states and their generated spectra of non-zero eigenval- 
ues, but also demonstrate that coherent optical informa- 
tion of multi-qubit can be stored in a medium with collec- 
tive effect in EIT. The storing and reading-out processes 
are controlled by stimulated photon transfer between the 
classical and quantized light field. In this section we will 
consider the transient dynamic process for the interac- 
tion between light fields and the EIT medium. Some ex- 
otic natures of light propagation in such a medium will 
be investigated with the help of the dynamic symmetry 
analysis. 

The evolution of the Heisenberg operators correspond- 
ing to the optical field and the collective excitations can 
be described as 

a k = -i\fNg k A k , 

Ak = -iVNg k a k - i£lC k , (36) 
Ck = -i£lA k 

for each light mode k. The above motion equations of 
a k , A k and C k with respect to a given optical mode k 
are coupled together, but there is no coupling of mode k 
operators to those of different optical mode k'. Due to 
the semi-direct product property of the dynamic group 

G = SU(2) <g> S <g> T, the generators T+, T_, and T 3 
of SU(2) algebra do not occur in the above system of 
Heisenberg equations. Here T is the Heisenberg- Weyl 
group generated by the creation and annihilation opera- 
tors a\ and a k . It leads to the evolution equations of the 
bright- and dark- state polariton operators 

Ak = -i^/g 2 k N + Q 2 B k , 

Bk = -iJg~lN + Q?A k , (37) 

D k =0 

in a straightforward way. The above equations mani- 
fest the basic features of the dark states: decoupled from 
other states and stable in the time evolution. 

With the initial conditions determined by those for C k 
and a k , the exact solutions of B k and D k are obtained 
as follows: 

D fc (t) = D fc (O) = a fc (O)cos0 fc -C fc (O)sin0 fc) (38) 
B fc (i) = O ie - lte « + 2 e it@k }, (39) 

01 = ^(B fc (0)-A fc (0)) 

2 = i(B fe (0) + A fe (0)) 



where a k (0), A k (0), B k (0) and C k (Q) are the initial 
Heisenberg operators, and Q k = \J 'g k N + Q 2 is the light 
field dressed Rabi frequency that measures the effective 
coupling of excitonic states to the external field. 

The dark-state polariton operator D k (t) is indeed a 
time-independent constant. That is, in the Heisenberg 
picture, the operator D k {i) always retains the initial 
k— mode dark-state polariton operator D k (0) for each 
k and the atomic system always keeps being an EIT 
medium. Here for simplicity we assume that the quan- 
tized light field propagates along the x— axis. Then, the 
positive frequency part of the quantized light field 

E + (x,t)=^J^a k (t)e ik *- ik * (40) 
k 

can be expressed explicitly in terms of the mode opera- 
tors of quantized light field 

a k (t) = D k (t) cos 9 k + B k (t) sin 6 k 

= [a k (0)cos9 k - C k (0)sm8 k ]cos8 k (41) 
-(O ie - 4 * 6fc +0 2 e ltek )sme k , 

which results from the expressions of bright- and dark- 
state polariton operators given above. 

With the above simple solution we can straightfor- 
wardly investigate how the quantized light propagates 
in the EIT medium. We assume that the collective state 
of the medium is initially in the ground state |b) and 
the initial state of the light field is a wave packet. The 
electromagnetic field quantized, the wave packet of light 
can be depicted by the direct product of many coherent 
states, or the multi-mode coherent state 

\a) =H®\a k ), (42) 
k 

where a k = exp[— f 2 (k — k ) 2 ] (up to a normalized fac- 
tor) since the probe light is a Gaussian wave packet in 
frequency domain. This is because the expectation value 
of the free quantized electromagnetic field operator 

i?+ (^ i ) = E v W afc(<)e " ct 

k 

resembles a classical wave packet as the superposition of 
infinite components of different frequencies. In interac- 
tion with EIT medium, the mean of a k (t) over the initial 
state |V>(0)) = |b) <g> \a) is 

(a k (t)) = (a fe (0))cos 2 6 k + 

sin^((O 2 )e- i * e '= + (O 1 )e it0fe ) 
= (a k (0)) (cos 2 9 k + sin 2 9 k cos Q k t) , (43) 

where we have used 



G 



(A fc (0)) = (C fc (0)) =0, 
(Bfc(O)) = (a fe (O))sin0 fe . 

So the mean of E + {x,t) is decomposed into three parts: 



(£+ M ) = £jM<a fc (0)> 



2eV 



(cos 2 6» fe + sin 2 9 k cosQ k t)e t{kx - kct) 
E+(x,t)+E+(x,t)+E+(x,t) (44) 



where 



Ei(x,t) = ^J^(a k (0))^9 k e^- 



a)-e k t] 



2eV 



E+(x,t) = ]T 2 J^(a k (0)) cos 2 fe e 



2 fl, -ife(x-ct) 



2eF 



(45) 



^>*) = E^ (°)> Sin2 B k e^- Ct ^ . 



2eV 



It should be noticed that the light field dressed Rabi fre- 
quency <d k modifies the dispersion relations of light prop- 
agation in the exciton medium. 

The three parts in the above decomposition of 
(E + (x,t)) can be understood as three wave packets 
spreading in the coordinate space. We can see this point 
by considering the initial quantized light field in the fre- 
quency domain as a Gaussian wave packet with the center 
frequency uj (= k n c = 0J a b)- For the first wave packet 
Ef(x,t), its position center is determined by maximizing 
the wave- vector-dependent phase: 



d 

— {kct + e k t-kx)\ ko =0. 



(46) 



This determines the group velocity of the wave packet 
E^(x, t) as 



9+ 



dx 



dk 



\k - 



(47) 



The term on the r.h.s. simply leads to a modification of 
the group velocity 



V g+ = c(l 



2/cqc y/l + n 



), 



(48) 



where n = g ^ 5 . This is one of the central results of this 
paper. It is evident that the first wave packet E+(x,t) 
represents the " super- luminal" light propagation with 
the group velocity larger than the light speed c. It is 
also emphasized that the " super- luminal" light propaga- 
tion can be understood according to the classical theory 
of wave propagation in an anomalous dispersion medium. 
It is the interference between different frequency compo- 
nents that results in this rather counterintuitive effect. In 



this sense, we think that the " super-luminal" light pulse 
propagation observed in experiments is not too odd and 
there is no direct connection between this phenomenon 
and the causality in relativity. 

It is not surprising to see the usual group velocity 
V go = c for the second wave packet Eo~(x,t). But for 
the third wave packet E^(x,t) we find the sub-luminal 
group velocity similarly: 



V„ 



dk 



= C(1 



2k c y/l + n 



)■ 



(49) 



This sub-luminal group velocity phenomenon is expected 
to find applications in quantum information processing. 
In fact, it is common sense to believe that the system 
stopping and slowing light could be used to store quan- 
tum information of photon qubits and in the storage 
time quantum information processing may be possible 
if we can sufficiently reduce the dissipative loss during 
the "reading" and "writing" operations. 

We can also analytically integrate out £^(x,i), 
E+(x,t) and E^(x,t) and obtain an explicit depiction 

of the propagating wave packet. Since g k = — p\J > 

it is convenient to write g k N as G 2 k with G = P\J 2 heV ■ 
Replacing the summation over k by integral and using 
(afc(0)) = a k , we obtain the normal part 



^0 



2L, 



ck n 2 



2eV n 2 + G 2 k 

xe -/ 2 (fc-fc )yfc(*-ct) dfcj ( 50 ) 



and the super-luminal and sub-luminal (or negative) 
parts 



E 



(+) 



ck G 2 k 



2eV n 2 + G 2 k 



where L me is the scale of integration along the propaga- 
tion direction (cc-axis) of the probe light. We then expand 

/ — O 2 'h 

the terms y/k, n -i_l G -± k , n^+c^fc an< ^ ® fe aroun d the cen- 
tral value ko of the wave vector of the input probe light. 
We can neglect the higher-order terms of A = k — ko since 
k is very close to ko. Then, with the defined parameters 
A = G 2 /£l 2 , Oo = + Ak 0} the approximate expres- 
sions of field components are obtained analytically as 



E ( +\x,t) 



and 



2/3(1 + Ako) 



{2f 2 + iDo[x-(c±Eo)t}} 



xexph^z- (c±E )t] 2 ]e ik °( x - c ±V (52) 
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-[2f 2 + iF (x - ct)] 



2/ 3 (l + Ak ) 1 

xexp[--^(x-ct) 2 ]e tko( - x - ct \ (53) 



4/2 



where c± = c ± ^ and 



Akn 



2k (l + Ak )- 

n A 
2{\ + Ak y 

1 -Ak 
2fc (l + Afc )' 



The above analytic result proves again the coexistence 
phenomena of both the super-luminal and sub-luminal 
(or negative) group velocities. From the Eqs. (52,53), 
it's obvious that the group velocity of the part Eq + \x, t) 
is V go = c, and the group velocity of wave packet E±(x, t) 
is 



V g± =c±E Q =c(l± 



2fc c y/T 



=) 



(54) 



which is the same as the Eqs. (48) and (49). 

Through the above theoretical analysis based on dy- 
namic algebraic method, several interesting properties of 
the light pulse propagation in an EIT medium are ob- 
served. The central result is the possibility of coexis- 
tence of "sub-luminal" and "super-luminal" group veloc- 
ity phenomena, occurring as a triple split of wave packet 
when spreading in such an EIT medium. We also observe 
that if we inject coherent light into the EIT medium for 
a certain time, there may appear three light pulses. One 
of these three light pulses still propagates at the veloc- 
ity c, but the other two have "super-luminal" and "sub- 
luminal" group velocities respectively. The " sub-luminal" 
group velocity might even have negative value in some 
case. This result is different from that of stable process in 
which only one single group velocity appears. This is be- 
cause we here only consider the transient phenomenon or 
property. It is also remarked that both "super-luminal" 
and "sub-luminal" (even or negative) group velocities are 
natural properties of light [2] and could be realized ex- 
perimentally. 



V. NUMERICAL SIMULATION OF WAVE 
PACKET SPLIT 

In this section we numerically simulate the dynamic 
process of the light pulse propagating in the above de- 
scribed EIT medium. According to the Eqs. (52,53), 
we calculate the shape evolutions of the wave packets 
of the three parts of the light field, respectively. Let 
^° = Ito = ^7 ^ e the central wavelength of the initial 



wave packet of the light pulse. To simulate the light pulse 
propagating in the EIT medium numerically we distin- 
guish the following two situations: (I) the spatial width d 
of light pulse is less than Ao; (II) d contains a few Ao's. 
In the second situation several oscillations will be clearly 
observed in a light pulse with a few frequency compo- 
nents. Let Ei(x, t) = E+(x, t)+ E^(x, t), for i = 0,+, — . 
For situation (I) that the spatial width d of light pulse 
is less than its central wavelength (d sa 0.3Ao), three 3- 
dimension curves of Eo(x,t), E+(x,t) and E-{x,i) are 
plotted as in Fig. 2. It is observed from Fig. 2 the three 
parts Eo(x,t), E + (x,t) and E_(x,t) of light field indeed 
have different group velocities. It's noted that in order 
to satisfy the near-resonance condition, the experimental 
probe light pluse contains large numbers of wavelengths 
much more than that given above where it's convenient 
to see the pulse splitting. 



(a) 




(b) 




(c) 




FIG. 2. Wave packet split of light pulse in the case that the 
spatial width of light pulse is less than its central wavelength 
d R3 0.3Ao): (a): the normal part with group velocity c. (b): 
the super-luminal part with group velocity being about 1.75c. 
(c) the sub-luminal part with group velocity being about0.25c. 

Fig. 3 is drawn for the situation (II) that the spatial 
width d of light pulse is large than its central wavelength ( 
d « 4Ao). In both of the situations it is observed that the 
group velocity (v ~ 1.75c) of E+(x,t) is larger than that 
of Eq(x, i), but the group velocity (v ~ 0.25c) of E-(x, t) 
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is less than that of Eo(x,i).The only difference between 
these two situations is the details of the oscillations of 
the caves. 



(a) 




FIG. 3. Wave packet split of light pulse in the case that the 
spatial width of light pulse is large than its central wavelength 
d pa 4Ao): (a): the normal part with group velocity c. (b):the 
super- luminal part with group velocity v — 0.25c. (c): a 
sub-luminal part with group velocity v = 1.75c. The other 
parameters is as the same as that in Fig. 2. 

When the system parameters are changed so that 



2k c y/l + n 

a seem-to-be exotic phenomenon of light pulse propa- 
gation is illustrated in Fig. 4. Here, while two parts 
of the splitting wave packet (Fig.4a-4b) possess nor- 
mal group velocity c and " super-luminal" group velocity 
(v ~ 2.36c) respectively, the third part propagates at a 
negative group velocity (v ~ —0.36c). 



(a) 




FIG. 4. Propagation of light pulse with a part possessing a 
negative group velocity (~ —0.36c). Here, the spatial width of 
light pulse is large than its central wavelength d « 4Ao): (a): 
the normal part with group velocity c. (b) the super-luminal 
part with group velocity v ~ —0.36c. (c) the part propagating 
with a negative group velocity (v ~ 2.36c). 

In fact the amplitude of these held components de- 
pends on the system parameters. It is easily seen from 
Eq. (45) that the intensity of E+(x,t) is equal to that 
of E-(x,t), but is not equal to that of Eo(x,t). The 
ratio of the amplitude of E±(x,t) to that of Eo(x,t) is 
approximately 

sin2 ^o _ 9jN _ n 
2cos 2 fco 2D 2 2' V ; 

When the ratio ^ is very small (<?; 1), most part of the 
light pulse will propagate in the EIT medium at the group 
velocity c. If the ratio ^ is very large 1), then most 
part of the light pulse propagates at " super-luminal" or 
sub-luminal (even negative) group velocity. To see the 
dynamic details clearly we also plot the caves of the spa- 
tial wave packets for different fixed instances in Fig. 5. 
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FIG. 5. Evolution of three splitting wave packets of light 
pulse. Here, d = 4A ) (a) The initial un-split wave packet with 
the spatial width of light pulse d « 4Ao);(b) [(c)] The shapes 
of the three splitting wave packets at times t\ ~ 2.3 * 10~ 13 s 
[ t 2 « 4.6* 10~ 13 s]. The cases with normal, "super- luminal" 
and sub-luminal are depicted respectively by solid line, dot line 
and dash line. The parameters are the same as that in Fig. 3 

Fig. 5 (a-c) shows the time evolution of these three 
wave-packets at t = to(= 0),ti,t 2 (0 < t\ < t 2 ) accord- 
ing to the Fig. 3. The field components E + (x,t) and 
E-(x,t)) with "super-luminal" and "sub-luminal" group 
velocities change their shapes of wave-packet during the 
evolution, but the part E (x,t) always propagates at c 
with an unchanged shape of wave-packet. Due to the 
intrinsic quantum coherence the EIT medium does not 
change the height of each wave-packet. Strictly speak- 
ing, c = c /no is not the vacuum velocity of light, while 
Co is,where no is the index of refraction [20] (here for 
simplicity, we set no = 1 in our numerical simulation in 
this work). To see the exotic phenomenon with negative 
group velocity, we take the 2-dimension curves at certain 
time in Fig. 6 according to the Fig. 4. Here, a inverse- 
direction "light propagation" can be seen clearly. 



FIG. 6. Phenomenon with a negative group velocity 
(v ~ —0.36c). Here, d = 4Ao): (a) The initial un-split 
wave packet with the spatial width of light pulse d w 4Ao);(b) 
[(c)] The shapes of the three splitting wave packets at times 
ti w 1.3 * l(n 13 s [ t 2 w 2.6 * lCr 13 s]. The cases with nor- 
mal and super-luminal are depicted respectively by solid line 
and dot line. The dash line describe the a part of the split 
wave packet propagating with a a negative group velocity 
(v ~ —0.36c). The parameters are the same as that in Fig. 4. 

We can also calculate the optical intensity in the case, 
where the initial quantized light field is multi-mode co- 
herent state. We have 

I(x,t) = (E-(x,t)E+(x,t)) 

= Eq Eq + E^E + + E^E_ 

+2Rc(E+E- + E+EZ + EZE~). (57) 

From the Eqs. (52) and (53), it's obvious that EqEq, 
E+ El and E+EZ are three approximate Gaussian wave 
packets respectively and E + EZ (also E+EZ and EZEq) 
is the interference term. At the initial time t — 0, I(x, t) 
is a Gaussian wave packet approximately since the center 
of each term contributing to I(x, t) is at x = 0. But after 
a certain time, the terms E + Eq , E+EZ and EZEZ arc 
nearly separated in spatial coordinate and the contribu- 
tion of interference terms will be close to zero. So in this 
case the intensity will only contain 3 separated Gaussian 
wave packets with the "super-luminal", "sub- (even neg- 
ative) luminal" and normal group velocities respectively. 

Finally we consider the quantum fluctuation of the 
probe quantum field. In the above discussions we assume 
the quantized light field of each mode to be initially pre- 
pared in its coherent state. This is too conceptual a setup 
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for experiment. If the initial state is a number state the 
expectation of the quantized light field operator vanishes 
for the representation of photon number conservation. 
For this reason we need to consider the quantum fluctu- 
ation described by the correlation function 

(E-(x,t)E+(x,t + T)) 
= £ hcl ^ )k e^(cos 2 9 k + sin 2 9 k cos Q k tf. (58) 

k 

It will determines the intensity spectra 

/oo 
dre-^ r (E-(x,t)E+(x,t + r)) = S(0,u). 
-OO 

(59) 

The one-time correlation is substantially the intensity of 
the probe light field 

(/(*,*)) =(E-(x,t)E+(x,t)) 

= J2 hC ^eV^ (C0S2 dk + Sln2 ° k C ° S Qkt)2 - m 

k 

As a modulated asymptotic function, this result is illus- 
trated in Fig. 7 for an initial number state with (n) k 
satisfying a Gaussian wave distribution in frequency do- 
main. It shows the instantaneous process of light prop- 
agation to approach a stable state. For very large Rabi 
coupling constant 0, 8 k ~ 0, Q k = J g k N + £1 2 ~ Q 
and the intensity approaches a constant (I(x, t)} — 
Efc hc ^}y )k . On the other hand for a very large medium 
enhanced coupling of probe light g k N or SI = 0, 9 k ~ 
f ,we can analytically show 



h c\k\ (n) k 



2eV 

^hc\k\(n) kc ^ 2{ 



(cos g 2 k Ntf 



2eV 



cp 2 N 
2freV 



kt), 



(61) 



which is asymptotic to a constant with the modulation 
frequency ui M = ^yk. 



1.5 



0.5 



FIG. 7. The intensity evolution of probe light with a mod- 
ulated asymptotic behavior. 




VI. CONCLUSION WITH REMARKS 

In conclusion, we have investigated how the quantized 
light propagates in a three-level A— type EIT medium 
by directly solving the Hcisenbcrg evolution equation of 
this light field based on the dynamic group method. For 
simplicity, the quantized light field is considered as res- 
onantly coupling to the A— type subsystem even though 
this light field contains a series of light modes and can 
not be resonant simultaneously. The physical reason 
for this consideration is that the light wave packet has a 
small width in frequency domain. In fact, the quantized 
light is considered as a quasi-plane wave and propagat- 
ing in the medium without the boundary effect. It should 
be noticed that the present treatment is valid only when 
the EIT medium is prepared in the low density excitation 
situation. 

We wish to emphasize again that it is also owing to 
the wave nature of light that the coexistence phenom- 
ena of both the "super- luminal" and " sub- luminal" (or 
negative) group velocities appears as predicted in this 
paper. It is very interesting to observe this coexistence 
phenomena experimentally and explore its potential ap- 
plication in quantum memory and quantum information 
process. To this end we need more details of physical 
considerations on the experimental techniques. For in- 
stance, we need to compare the size of the sample of the 
EIT medium and the split distances of three wave pack- 
ets resulting from the light pulse. In principle due to the 
different group velocities the evolution of sufficiently-long 
time will distinguish the wave packets, which might not 
preserve their shape because of wave packet spreading or 
the dissipation and decoherence due to the coupling to 
environment. On the other hand, what we predict are 
only the transient phenomena. So it is somehow diffi- 
cult to observe this coexistence phenomena of both the 
" super- luminal" and " sub-luminal" (or negative) group 
velocities in a practical experiment. 
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APPENDIX: LIGHT PROPAGATION IN 
EXCITON SYSTEM: ANALYTIC RESULT FROM 
SUB-DYNAMICS MEHTOD 



In general we consider a physical system with a dy- 
namic symmetry characterized by a Lie group G. This 



means that the Hamiltonian of the considered system 

H = H[G] = H(g u g 2 ,...) (A.l) 

is a functional of the generators g\, g 2l --- of G. These 
generators can be understood as the basic dynamic vari- 
ables of the system. Suppose there exists a subgroup 
S C G such that 

[H, G] C S. (A.2) 

Let si, S2, ... be the generators of S. Then the system of 
the Heisenberg equations 

i^s k = [s k ,H]e S (A.3) 

about si, s 2 ,... are closed to form the so-called sub- 
dynamics. Through the subset {si, s 2 , ...} of the com- 
plete set of dynamic variables gi, g 2 , the sub-dynamics 
depict the main features of the considered system. 

The excitonic system in this paper serves as a practical 
example for 

G = SU(2)®~®r (A.4) 

and S = S ® T where T is the Heisenber-Weyl group 
generated by the creation and annihilation operators a' k 
and a/j of light field. 
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